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Q_i' Color-octet contributions to the associate J/'i/' + 7 hadroproduction are stud- 
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ied in detail, and found to be negligible, compared to the ordinary color-singlet 
contribution. Within the color-singlet model the J/V' + 7 production in the 



^ ■ leading order is possible only through gluon-gluon fusion process. Therefore, 



the associate J/ip -\- 'j hadroproduction remains to be useful as a clean chan- 
nel to probe the gluon distribution inside proton, to study heavy quarkonia 
production mechanism, and to find proton's spin structure. 
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I. INTRODUCTION 



Up to very recently the conventional color-singlet model had been used as only pos- 
sible mechanism to describe the production and decay of the heavy quarkonium such as J/ ijj 
and T. An artificial K-factor (~ 2-5) was first introduced to compensate the gap between 
the theoretical prediction and experimental data. The uncertainties related to this large 
i^T-factor are the following: We don't know precisely the correct normalization of the bound 
state wave function, the possible next-to-leading order contributions, the mass of the heavy 
quark inside the bound state, etc. However, even with this large K-factor some experimental 
data are found to be difficult to describe. For the direct J/ ip production with large in pp 
collisions, the dominant mechanism has been found to be through final parton fragmentation 
. As applications of this fragmentation mechanism, various studies of prompt charmonium 
production at the Tevatron collider have been carried out, and the CDF data on the 
prompt Jjip production qualitatively meet these theoretical predictions. Nevertheless, 
the ip' production rate at CDF is about 30 times larger than the theoretical predictions, 
which is the so-called %p' anomaly, even after considering the fragmentation mechanisms of 
g ip' and c ip'. As a scenario to resolve this ip' anomaly, the color-octet produc- 
tion mechanism was proposed |^. By using this idea, heavy quarkonium (charmonium) 
hadroproduction through the color-octet (cc)8 pair in various partial wave states {^^^^Lj) 
has been studied in addition to the color-octet gluon fragmentation approach HJl^l- We 
note that the inclusive T production at the Tevatron also shows an excess of the data over 



theoretical estimates based on perturbative QCD and the color-singlet model |Tl|. In this 
case, the Pj, of the T is not so high, so that the gluon fragmentation picture may not be a 
good approximation any more. 

Since it has been proposed that the color-octet mechanism can resolve the ^' anomaly 
at the Tevatron collider, it is quite important to test this mechanism at other high energy 
heavy quarkonium production processes. Up to now, the following processes have been 
theoretically considered: inclusive J/ip production at the Tevatron collider and at fixed 
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target experiments [§,0,0], spin alignment of leptons to the decayed J /if) [|T3[, the polar 
angle distribution of the J /if) in e+e^ J /ip + X inclusive J /if) productions in B meson 
decays [|15],|l6| and in decays at LEP |T^-|1^], J/ip photoproduction p^ , |20| -p^, and color- 



octet J/ip production in T decays ^3[. For more details on recent progress, see Ref. p3| , 
which summarizes the theoretical developments on the quarkonium production. Recently, 
the helicity decomposition method in NRQCD factorization formalism was developed by 
Braaten and Chen With this method, polarized J/ip production in B decay was 

considered in Ref. In Ref. pSj^ , it was pointed out that there are interference terms 



of different ^Pj contributions in polarized J/ip production. 

In this work, we investigate the color-octet mechanism for associate J /il) + production 
in hadronic collisions. The associate J /ip + production has been first proposed as a clean 
channel to probe the gluon distribution inside the proton or photon [0], and then to study 



the heavy quarkonia production mechanism [29| as well as to investigate the proton's spin 



structure [0]. If the color-octet mechanism gives a significant contribution to this process, 
the merit of this process to probe the gluon distribution etc. would be decreased or one 
should find suitable cuts to get rid of this new contribution. Associate J /ip + production 
has also been studied as a significant QCD background to the decay of heavier P-wave 



charmonia {xcj{P)) in fixed target experiments ||3T1 , |32| within the color-singlet model ^3 
In Ref. 0], the fragmentation contribution from p + p {coTg)+'~f + X^J/tp + 'j + X 
at Tevatron energies {^/s=1.8 TeV) is studied, and these fragmentation channels are found 
to be significantly suppressed compared to the conventional color-singlet gluon-gluon fusion 
process. 

In Section 2, we explain in detail the method to deal with inclusive heavy quarkonium 
production, following the nonrelativistic-QCD (NRQCD) factorization formalism given by 
Bodwin, Braaten and Lepage (BBL) [^]. A summary of the kinematics related to the p^ and 
Mj/^^^ distributions is also given in this Section. In Section 3, we discuss in detail the p^ 
and Mj/^+^ distributions in the hadronic J/ip + ■y production, and compare the color-octet 
contribution with the color-singlet gluon-gluon fusion contribution. Section 3 also contains 



our conclusions. 



II. HADROPRODUCTION OF J/V' + 7 

A. General Discussions 

Now we consider the associate production oi & J/ip and a photon at hadronic colhders 
and fixed target experiments, including the new contribution from the color-octet mechanism 
in addition to the usual color-singlet contribution. The only possible subprocess in the 
framework of the color-singlet model in the leading order is through gluon-gluon fusion, 

9 + 9-^1 + (cc)(35{'))(^ J/^). (See Fig. 1) (2.1) 

This gluon-gluon fusion process for the associate 7/?/^ + 7 production is again possible within 
the color-octet mechanism. (See Fig. |I]). However, in the leading order color-octet contri- 
butions there also exist many more subprocesses, which are 

W), (2.2) 

J/^), (2.3) 

W), (2.4) 

J/^), (2.5) 

J/^). (2.6) 

These subprocesses are possible through the effective vertices of 

^ + g^{cc){^SP or =^Pf), (2.7) 

which is shown in Fig. |], and 

g + g-(cc)(^5f)), (2.8) 

which is shown in Fig. ^. We note that the (^S'j^'^'') channel in Eq. ( |2.7| ) is vanishing, even 
though the gluon is not on its mass shell. The Feynman diagram for the gluon initiated 
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g + g^^+{cc){'Sf^){- 

g + g^7+(cc)(^4'))(- 
g + g-7+(cc)(M'^)( 
g + g-.7+(cc)(^5f))(- 



subprocesses, Eqs. (p.2| , |27^ ), is shown in Fig. ^ (a). For the case of quark initiated 
subprocesses, there are two kinds of diagrams; The Feynman diagram for quark initiated 
subprocesses Eqs. ( p.4| , p.5| ), due to the effective vertex Eq. (p. 7]), is shown in Fig. ^ (b). 



The other diagrams, which are possible through the effective vertex, Eq. ( |2.8| ), are shown 
in Fig. |[ 

We first perform a naive power counting analysis of those various subprocesses. Free 
particle amplitude-squared for the color-singlet gluon-gluon fusion process is [C(aaf)]. 
In the NRQCD framework the color-octet (QQ) states can also form a physical J/ijj 



state with dynamical gluons inside the quarkonium with wavelengths much larger than the 
characteristic size of the bound state {^^1 / (MqVq)) . In Coulomb gauge, which is a natural 
gauge for analyzing heavy quarkonium, these dynamical gluons enter into the Fock state 
decomposition of physical state of J/ip as 

I JM = OmQQ)es[''>)) + 0{vq)\{QQ)Cp!j%) 

+ 0{vl)\{QQ)Cs[''%g) + 0{v'QmQ)Cs!,%) + ... . (2.9) 

In general, a state {QQ)(^^^^Lj) can make a transition to {QQ)(^^^^{L ± l)j), or more 
specifically {QQ)(^Pf) {QQ){^S[^^) through the emission of a soft gluon (chromo- 
electric dipole transition), and it is an order of Vq suppressed compared to the color-singlet 
hadronization. For the case of chromo-magnetic dipole transition, such as {QQ){ Sq ) — >■ 
{QQ)(^S[^^), it is suppressed by vq at the amplitude level. If we also consider the fact that 
the P wave state is MqVq order higher than S wave state in the amplitude level, one can 
naively estimate that the transitions 

{QQ){^Pf^)^ J/ij + X and (QQ)(i5f ) ^ J/^ + X (2.10) 

are commonly Vq order suppressed at the amplitude-squared level, compared to the transi- 
tion 

(gg)(35«)^ j/^. (2.11) 
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Since these color-octet processes have the same order 0{aa'^) as that of the color-singlet 
gluon-gluon fusion process in free particle scattering amplitude, the color-octet subprocesses 
are order Vq suppressed compared to the color-singlet gluon-gluon fusion subprocess. 

However, such analyses are only naive power counting, and do not guarantee that the 
color-octet contributions are suppressed compared to the color-singlet gluon-gluon fusion 
contribution all over the allowed kinematical region. If we consider inclusive J/ip photo- 
production via 2^2 subprocesses, as shown in Refs. |]16|,^, color-octet contributions 



dominate in some kinematical region, even if the naive power counting predicts the sup- 
pression of color-octet contributions compared to the color-singlet one. In this respect, it is 
worthwhile to investigate in detail how much the color-octet mechanism contributes to the 
J/V' + 7 hadropro duct ion. 

In order to predict numerically the physical production rate, we need to know a few 
nonperturbative parameters characterizing the fragmentation of the color-octet objects into 
the physical color-singlet J/ip. These nonperturbative matrix elements for the color-octet 
operators have not been determined completely yet. After fitting the inclusive J/ip pro- 
duction at the Tevatron collider, using the usual color-singlet J/ production, the cascades 
production from xd^P), and the new color-octet contributions, the authors of Ref. |]TU[ have 
determined 

(0|O^(^5i)|0) = (6.6 ± 2.1) X 10"^ GeV^ (2.12) 



M2 3 

with Mc = 1.48 GeV. Although the numerical values of two matrix elements, {0\Og{^Po)\0) 
and {0\Of Sq)\0) , are not separately known in Eq. (|2.13|) , one can still extract some useful 
information from them. Assuming both of the color-octet matrix elements in Eq. ( p.l3| ) are 
positive definite, then one has [] 



^The matrix element {0\O^ Po)\0) can be negative whereas (0|C'g("^ 5*0)10) is always positive 
definite [26|. Here we choose these ranges just for simplicity. 
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< {0\OtCSo)\0) < (6.6 ± 1.5) X 10-^ GeV^ (2.14) 
< (^\^s^P<^m ^ (2.2 ± 0.5) X 10"2 GeV^ (2.15) 

These inequalities could provide us with a few predictions on various quantities related to 
inclusive J/ip productions in other processes, and also enable us to test the idea of color-octet 
mechanism in the associate J/ip + ^ production process. 



B. Kinematics 

Let us consider the process 

«(Pi)mpi) + b{P2)/p{P,) ^ J/i^{P) + l{k). (2.16) 

in pp (or alternatively pp) collisions. We can express the momenta of the incident hadrons 
(p, p) and partons (a, h) in the pp CM frame as 

Pi = ^(1, +1, 0), Pi = ^(^1, +^1, 0), (2.17) 



^2 = ^(1, -1, 0), P2 = ^(2^2, -X2, 0), (2.18) 

where the first component is the energy, the second is longitudinal momentum, and the third 
is the transverse component of the particle's momentum. The variables X\ and ^2 are the 
momentum fractions of the partons. The momenta of the outgoing particles are given by, 

P = {E^, Pi, +Pt) = (Mt cosh y^, Mrsinhy^, +Pt), (2.19) 
k = {E\PI-Pt) = ( PTCOshy^ PTsinhy^ -Pt), (2.20) 

where Pr is the common transverse momentum of the outgoing particles, Mr is the trans- 
verse mass of the outgoing J/il^, and y'^ (or y'^) is the rapidity of J/V' (or 7). Some useful 
relations among the above variables are: 

E"^ + Pt = MTe+y\ E^-Pt = MTe-y\ (2.21) 
E^ + P2^ Pre+y", E^-Pl^ Pre-y". (2.22) 



The Mandelstam variables are defined respectively as 

s = (Pi + P2)', (2.23) 

s = (pi + = 2X1X2P1 • P2 = (2.24) 

£ = (pi - A;)2 = -x^^sPre-y' = {p^ - Pf ^ - x^^sMTe^^^ (2.25) 

u^{p2- = -X2A/iPTe+^" = (pi - Py = M| - xiv^Mre-f*, (2.26) 

M| = s + t + 'a. (2.27) 

The last equation, representing the energy-momentum conservation, leads to 

v^MT(x2e^* + xie-^*) = s + M|, 

v^Pt (a^se^' + xie-^') = s - M}. (2.28) 
After introducing dimensionless variables, 

xt = 2Pt/Vs, xt^2Mt/Vs and r = M^/s, (2.29) 

we can solve the energy momentum relation to get X2 in terms of other variables as 

XiXtC^^^ - 2t xixre-y'' + 2r 

2xi-x,e+^- = 2x, - xre^y-^ " ^^'^''^ 

We also obtain the relations for rapidities of the outgoing particles in terms of Xi, X2, xt, 
xt and r as following 

{xiX2 + r) + ^J{x^x^'+Ty'^^^^X2^ 



X2XT 



, ^. {X1X2 -t)-^ {xiX2 - t)2 - XiX2xl 

exp = . 2.31) 

X2XT 

In order to get the distributions in the invariant mass Mj/^_^.y{= and transverse mo- 
mentum Pj, for the process g + g ^ J/ip + j process, we express the differential cross section 
as 

= fg/pi^i: Q^)fg/p{x2: Q'^)^dxidx2di 

= f9/p{xi,Q'^)fg/p{x2,Q'^)^J dxidxTdMj/^+^ 

= fg/p{xi,Q')f,dx2,Q')^J (^^] dx.dy^dPr, (2.32) 



where the corresponding Jacobians are given by 

j( \ _ 2x,xjMj/,^, j(^^]= (2.33) 

\xiXTMj/^+yJ XT{x2e^y* — xie y^) \xiy'^PT/ 2xi — xre+y 

Then the distributions are expressed as 

= / dx, dxrJ ( ] , (2.34) 

dMj/^+^ J \XiXtMj/^+^J dxidx2dt 

r,,^iy*jljilE^\_f^^, (2,35) 



dPr J \xiy'^ Pt ) dxidx2dt 

And the allowed regions of the variables are given by 

< < 1, 

s 

r2 „ ^ /\/r2 



< XT < ^ , . (2.36) 

^JXiX2 



C. The Nonrelativistic-QCD (NRQCD) Factorization Formalism 

First we consider the general method to get the NRQCD cross section for the process 
a + h ^ {QQ){^^^^Lj^'^^){—>- H) + c, where H is the final state heavy quarkonium and 
(QQ)(^^+^^j'^^) is the intermediate (QQ) pair which has the corresponding spectroscopic 
state. From now on, we use the subscript n to represent the spectroscopic (QQ) state of 
(25+i^{^"-i.8)^^ for simplicity. Once the on-shell scattering amplitude of the process A{a+b — *■ 
Q + Q + c) is given, we can expand the amplitude in terms of relative momentum q of the 
quarks inside the bound state because the quarks, which make up the bound state, are heavy. 
For more details of the method to deal with the heavy quarkonium production following the 



BBL formalism |3^, we refer to Refs. P, pl)| , p!6| . The heavy quarkonium H production cross 



section a{pi) + b{p2) {QQ)n{P) + c{p^) ^ H + c + X is given by 

5(a(pi) + b{p2) ^ {QQ)n{P) + c(p3) ^H + c + X 
dt ^ 



da' (OlOflO) 
X ^ X ^ ' " ' , (2.37) 



dt 2J 



where 

da' 1 



-^|A^' (a(pi) + b{p2) ^ (QQUP) + c{ps)) \'. (2.38) 



dt IGnP 
Here, M.' is the amphtude of the process 

a(pi) + b{P2) ^ (gg)„(P) + c(p3), (2.39) 

which can be obtained by integrating the free particle amphtude over the relative momentum 
of the quark inside the intermediate state {QQ)n{P), after projecting appropriate spectro- 
scopic state Clebsch-Gordon coefficients. The parameter Cn is defined by 

{2Nc (color — singlet), 
(2.40) 
- 1 (color - octet). 

And (0|(!?^|0) is the non-perturbative matrix element representing the transition 

{QQ)C'+'L^j^'^)^H. (2.41) 

Finally, J denotes the angular momentum of the intermediate state {QQ){^^^^L^j'^^), not 
of the physical state H. For the case of color-singlet intermediate state, which has the 
same spectroscopic configuration with H, we can relate the matrix elements to the radial 
wave-function of the bound state as 

'^|i?5(0)P (5 -wave), 

(2.42) 



Cn X (2 J + 1) 



|-|i?^(0)p (P-wave). 



For example, if we consider the ^jj production via (^,5^ ^), (^^f ^), {^S?^), (^Pq^^^), i^Pi^^) and 
(^P2^^) intermediate states, then the partonic subprocess cross sections are given by 

S(„,tet) = ^ ( ^('Sf ) . (0|O*CSf )|0> + £(»SP') X 

dt^ ' SM^ydt^ ° ' ^ ' ^ ^1 / dt^ ' 3 
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after imposing the heavy quark spin symmetry 

(0|C-^/^(=^Pf )|0) = (2 J + l)(0|0-^/^(=^Po^'^)|0). (2.44) 
The subprocess cross section for the color-singlet gluon-gluon fusion is well known: 



da . , Ml 

where the overall normalization Mi is defined as 

4 



s - 4:M^y{i - AM^f{u - 4M2)2 



(2.45) 



Mi = ^ (47ra,)^(47ra)e^ M^Gi{J/ij). (2-46) 
9 



The parameter Gi{J/ip), which is defined in NRQCD as 

GiiJm ^1^5(0)1 , (2.47) 

is proportional to the probability of a color-singlet (cc) pair in the {^S[^^) partial wave state 
to form a physical J/ ip state. It is related to the leptonic decay width 

r(J/V' ^ l^r) = ^ nela^ Gi(J/^), (2.48) 

where Cc = 2/3. From the measured leptonic decay rate of J/ip, one can extract 

GiiJ/ip) ^ 106 MeV. (2.49) 

After including the radiative corrections of 0{as) with as{Mc) = 0.27, this value is increased 
to ~ 184 MeV. Relativistic corrections tend to increase Gi{J/ip) further to ~ 195 MeV [jl5[. 



For the case of color-octet subprocesses, we can use the average- squared amplitude of 
the processes 7 + 5' (or g) {cc){^^^^L^f') + g (or q) in Ref. |]T6|, after crossing {k — »• 
-k and q2 -^2) 



l\Mfij + g^ {cc)C'^'Lf)+g){i,s,u), (2.50) 

25•+l^{8)^ 



Y^lMfiq + {cc)C'+'Lf) + ^){s, t, u) 



\\M'\\^ + q^ {cc)C'+'Lf) + q){i,s,u). (2.51) 
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As previously explained, there is only one color-singlet subprocess, but 9 color-octet 
subprocesses are contributing to the process 

P + V{P)^ Jl^ + l + X. (2.52) 

Whereas the initial partons for the color-singlet process are only the gluons, quarks can also 
be the initial partons for the color-octet subprocesses. For the (^S'j^'') channel, the gluon 
contribution corresponding to Fig. ^ (a) is absent because the effective vertex corresponding 
to Fig. ^ is vanishing. Using the parton level differential cross sections for various channels, 
in the next section we discuss in detail the and Mj/^+^ distributions in the hadronic 
J/?/' -|- 7 production. We also compare the color-octet contributions with the color-singlet 
gluon-gluon fusion contribution. 



III. NUMERICAL RESULTS AND CONCLUSIONS 

Now, we are ready to show the numerical results from the analytic expressions obtained 
in the previous section. The results are quite sensitive to the numerical values of QCD 
coupling constant a^, mass of charm quark M^, and the factorization scale Q. For numerical 
predictions, we used ^.(M^) = 0.27, = 1.48 GeV and = (Mr/2)2, where Mr is the 
transverse mass of the outgoing J/ip. For the structure functions, we have used the most 
recent ones, MRSA and CTEQ3M ||^, and it turns out that both sets of structure 



functions give more or less the same results within ~ 10%. 

We use the numerical value of the color-octet matrix element, shown in Eq. (2.12), 

{Q\O^Csf)\0) = (6.6 ± 2.1) X lO^^GeVl (3.1) 

And since the matrix elements {0\Ot So)\0) and {0\Ot Po)\0) are not determined sepa- 
rately, we present the two extreme values allowed by Eqs. (2,13,2. 14, |2.15| ) as 

^5*0 saturated case : {0\OtCSo)\0) = (6.6 ± 1.5) x 10"^ GeV^ 
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^Pj saturated case : (0|C^(3Po)|0) = (2.2 ± 0.5) x lO'^ GeV^ x M^, 

(0108^150)10) = 0. (3.2) 

After analyzing the numerical results, we found that the color-octet contribution of ('^S'l'^^) 
was significantly suppressed compared to the others over the entire phase space. Therefore, 
we do not present the result for the channel (^S'|^'*). 

We first consider the Mj/^+^ distribution at a fixed target experiment. In Ref. p3 



Berger and Sridhar investigated the J/'?/' + 7 invariant mass distribution at ISR and 
E705 [0 experiments. They argue that this process forms an important (and computable) 
part of the background to the decay of P-state charmonia into J/^/^ + 7 at low values of the 
invariant mass of the J pair. In Fig. ^, we present the Mjj^^^ distribution for the E705 
fixed target experiment (a/s = 23.72 GeV). We note that the value of Mj/^+^, to produce 
the J/?/' + 7 pair, must be larger than Mj/^. The solid line is the color-singlet gluon-gluon 
fusion contribution. The dashed line represents ^Pj^''-saturated color-octet contribution, and 
the dotted one is 5*0 -saturated color-octet contribution. As shown in this figure, the color- 
octet contributions are strongly suppressed compared to that of the color-singlet gluon-gluon 
fusion process. Considering the fact that the ^pj^''-saturated and ^S'o^'' -saturated results are 
maximally allowed (upper bound) contributions, as explained in Eqs. (2.14,2.15,3.2), the 
color-octet contributions must be smaller than the color-singlet one by more than one order 
of magnitude all over the region. 

In Ref. the color-singlet gluon-gluon fusion contribution has been compared with the 
gluon and charm-quark fragmentations to inclusive J/ production at Tevatron energies. It 
was found that the fragmentation contributions are more than an order of magnitude smaller 
than the color-singlet gluon fusion contribution over the whole range of considered. In Fig. 
1^, we show Bda jd^rj, for J /ip+'~^ production in pp collisions at -y/i =1.8 TeV, integrated over 
a rapidity interval \yj^^\ < 0.5, where B is the J/ip branching ratio into leptons (P=0.0594). 
All the line styles in Fig. 7 are identical with those of Fig. |^. As shown in Fig. 0, the 
color-singlet component is dominant where < Q GeV. There is a cross over point where 
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color-octet components dominate that of color-singlet gluon-gluon fusion subprocess around 
= 6 GeV. Though the color-octet contributions dominate in the high region, the cross 
section in this region is much smaller than that of low p^, region, where the color-singlet 
contribution dominates. As previously explained, the ^Pj^''-saturated and ^^o^^-saturated 
results are maximally allowed (upper bound) contributions, and the color-octet contributions 
might be smaller than the color-singlet one even at large p^, > 6 GeV. 

In Fig. Hwe show the rapidity distribution of J/ip, Bda/dy, integrated over the transverse 
momentum of J/ip, 3 GeV < p^, < 6 GeV. The rapidity distribution integrated over this 
region is dominated by the color singlet gluon fusion process, as shown in Fig. ^ Considering 
the fact that the matrix elements (0|C'^/'^(^5o)|0) and (0|0-^/'^(3Po)|0) in Eq. Q are pos- 
sibly overestimated by about an order of magnitute, the color octet contributions shown in 
Fig. H could be negligible compared to the color singlet gluon fusion process. Finally, we note 
that even though we have used the covariant projection method in our calculation, we have 
considered only the unpolarized J/ip production. Therefore, there are no interference terms 
among different ^Pf'^ components, and our results shown here are perfectly valid P5| , p7[ . 

To conclude, we have considered the color-octet contributions to the associate J/ip + 
7 production in the hadronic collisions, also compared to the conventional color-singlet 
gluon-gluon fusion contribution. Within the color-singlet model the J/ip + production 
in the leading order is possible only through gluon-gluon fusion process. As we expected 
according to the naive power counting, we found that the color-octet contributions are 
significantly suppressed compared to the color-singlet gluon-gluon fusion process. For the 
case of the hadroproduction at the Tevatron energies, we found that there is a cross-over 
point at high p^ where the color-octet contributions become dominant, and it is still larger 
than the fragmentation contributions shown in Ref. ||3^. This suppression of the color- 
octet contributions in the associated production is qualitatively consistent with the recent 
result of Cacciari and Kramer |^ . They investigated the J/?/) + 7 production via resolved 
photon collisions in HERA, and found the strong suppression on the color-octet contributions 
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compared to the color-singlet one as at 



100 GeV for pr > 1 GeV |22 




15% 



(3.3) 



Though most events at the Tevatron collider are in the region where the color-singlet gluon- 
gluon fusion contribution dominates, one could require an additional cut, for example: < 
6 GeV, to guarantee that the color-singlet gluon-gluon fusion process remains one of the 
cleanest channels to probe the gluon distribution inside proton, and to study heavy quarkonia 
production mechanism. 
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FIGURE CAPTIONS 



Fig.l Feynman diagrams for the color-singlet subprocess ior g + g ^ (cc){^S[^^) + 7 and 
the color-octet subprocess ioi g + g ^ {cc){^Sq^^ or ^Pj^"*) + 7. 

Fig. 2 Feynman diagrams for the effective ^ — g — {cc){^Sq^\ or ^Pjf^) vertex. 

Fig. 3 Feynman diagram for the effective q — q — {cc){^S[^^) vertex. 

Fig. 4 Feynman diagrams for the color-octet contributions to the subprocesses {&) g + g ^ 
{cc)CS^'^\ or ^Pf) +7 and (b) g + g ^ {cc)Csi^\ or ^pj^^) + 7. 

Fig. 5 Feynman diagrams for the color-octet contributions to the subprocess g + g — > 
(cc)(35f))+7. 

Fig. 6 Invariant mass Mj/^_|_^ distribution of the process pp J/tp + where y/s — 
23.72GeV. The sohd hne is the color-singlet gluon-gluon fusion contribution. The dashed 
hne represents the ^Pj^^-saturated color-octet contribution, and the dotted one is the ^Sq^^- 
saturated color-octet contribution. 

Fig. 7 Transverse momentum of J/ip, Pj,{J/ip), distribution Bda/dpj,, integrated over the 
J/t/j rapidity range \yj^^ \ < 0.5, for the process pp J/t/j+j+X, where i/i — 1.8 TeV. Here, 
B is the branching ratio of J/ijj decay into leptons (P=0.0594). The solid line represents the 
color-singlet gluon-gluon fusion contribution. The dashed line represents the ^Pj*^ -saturated 
color-octet contribution, and the dotted one is the ^S",]^^ -saturated color-octet contribution. 

Fig. 8 Rapidity of J/t/j, y{J/'4>), distribution Bda/dyj^^, integrated over the J/ip trans- 
verse momentum range 3 GeV < p^, < 6 GeV, for the process pp ^ J/ip + ^ + X , where 
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^/s — 1.8 TeV. Here, B is the branching ratio of J/V' decay into leptons (5=0.0594). The 
sohd hne represents the color-singlet gluon-gluon fusion contribution. The dashed line repre- 
sents the ^Pj^''-saturated color-octet contribution, and the dotted one is the ^5'o^'* -saturated 
color-octet contribution. 
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FIGURES 
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FIG. 6. 
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FIG. 7. 
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